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ABSTRACT 

We study the expected redshift evolution of galaxy cluster abundance between z ^ 3 in different 
cosmologies, including the effects of the cosmic equation of state parameter w = p/p. Using the halo 
mass function obtained in recent large scale numerical simulations, we model the expected cluster yields 
in a 12 deg^ Sunyaev-Zel'dovich Effect (SZE) survey and a deep 10'* deg^ X-ray survey over a wide 
range of cosmological parameters. We quantify the statistical differences among cosmologies using both 
the total number and redshift distribution of clusters. Provided that the local cluster abundance is 
known to a few percent accuracy, we find only mild degeneracies between w and either or h. As 
a result, both surveys will provide improved constraints on f2m and w. The ilm-w degeneracy from 
both surveys is complementary to those found either in studies of CMB anisotropics or of high-redshift 
Supernovae (SNe). As a result, combining these surveys together with either CMB or SNe studies 
can reduce the statistical uncertainty on both w and fi^ to levels below what could be obtained by 
combining only the latter two data sets. Our results indicate a formal statistical uncertainty of « 3% 
(68% confidence) on both and w when the SZE survey is combined with either the CMB or SN 
data; the large number of clusters in the X-ray survey further suppresses the degeneracy between w and 
both Qm and h. Systematics and internal evolution of cluster structure at the present pose uncertainties 
above these levels. We briefly discuss and quantify the relevant systematic errors. By focusing on 
clusters with measured temperatures in the X-ray survey, we reduce our sensitivity to systematics such 
as non-standard evolution of internal cluster structure. 
Subject headings: cosmology: theory - cosmology: observation 



1. INTRODUCTION 



It has long been realized that clusters of galaxies provide 
a uniquely useful probe of the fundamental cosmological 
i paramotorD. — The formation of the largo ooalo dark mat — 



tor (DM) potontial woUq of oluDtorp io likoly indopondont 
of complex gas djrnamical processes, star formation, and 
feedback, and involve only gravitational physics. As a re- 
sult, the abundance of clusters A^tot and their distribution 
in redshift dN/dz should be determined purely by the ge- 
ometry of the universe and the power spectrum of initial 
density fluctuations. Exploiting this relation, the observed 
abundance of nearby clusters has been used to constrain 
the amplitude dg of the power spectrum on cluster scales 
25% (c.g 



to an accuracv of - 



1993, Viana fc Liddle 1996) 



White, Efstathiou & Frenk 
'I'he value of in these 
studies depends on the assumed underlying cosmology, es- 
pecially on the density parameters flm. and ^Ia. Subse- 



quent works ( 


Bahcall & Fan 1998 




Blanchard & Bartlett 


1998 




Viana & Liddle 1999) have shown that the redshift- 



constrains on these two cosmological parameters. 

In the above studies, the equation of state for the A- 
component has been implicitly assumed to be p = wp 
with w = —1. The recent suggestion that w might be 
different from —1, or even redshift dependent (Turner & 



White 1997| , Caldwell, Dave fc Steinhardt 1998D has in- 
spired several studies of cosmologies with a component of 
dark energy. From a particle physics point of vie w, such 
w > ~1 can arise in a number of theories (see [Freese, 
Adams fc Frieman 1987|, |Ratra fc Peebles 1988i |l'urner & 
White 1997| , [Caldwell, Dave fc Steinhardt 1998| and refer- 
ences therein)^ It is therefore of considerable interest to 
search for possible astrophysical signatures of the equation 
of state, especially those that distinguish w = — 1 from 
w > —1. Wang et al. (2000) has summarized current as- 
trophysical constraints that suggest —1 < w < —0.2; while 
recent observations of Type la SNe suggest the stronge r 
constraint w ^ —0.6 ( Perlmutter, Turner fc White 1999 ). 

The galaxy cluster abundance provides a natural test 
of models that include a dark energy component with 
w ^ —1, because w directly affects the linear growth of 
fluctuations D^, as well as the cosmological volume ele- 
ment dV/dzdfl. Furthermore, because of the dependence 
of the angular diameter distance dA on w, the experimen- 
tal detection limits for individual clusters, e.g., from the 
Sunyaev-Zel'dovich effect (SZE) decrement or the X-ray 
luminosity, depend on w. Wang fc Steinhardt (1998, here- 
after WS98) studied the constraints on w from a combina- 
tion of measurements of the cluster abundance and Cos- 
mic Microwave Background (CMB) anisotropics. Their 
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work has shown that the slope of the comoving abun- 
dance dN/dz between < z < 1 depends sensitively on w, 
an effect that can break the degeneracies between w and 
combinations of other parameters {h, 51, n) in the CMB 
anisotropy alone. 

Here we consider in greater detail the constraints on 
w, and other cosmological parameters, from cluster abun- 
dance evolution. Our main goals are: (1) to quantify the 
statistical accuracy to which w ^ —1 models can be distin- 
guished from standard A Cold Dark Matter (CDM) cos- 
mologies using cluster abundance evolution; (2) to assess 
these accur acies in two specific c luster surveys: a deep 
SZE survey ( Carlstrom ct al. 1999D and a large solid angle 
X~ray survey, and (3) to contrast constraints from cluster 
abundance to those from CMB anisotropy measurements 
and from luminosity distances to high-redshif t Supernovae 
dSchmidt ct al. 1998| , [Pcrlmuttcr et al. 1999| ). 

Our work differs from the analysis of WS98 in sev- 
eral ways. We examine the surface density of clusters 
dN/dzd^l, rather than the comoving number density n{z). 
This is important from an observational point of view, be- 
cause the former, directly measurable quantity inevitably 
includes the additional cosmology-dependence from the 
volume element dV/dzdQ. We incorporate the cosmology- 
dependent mass-limits expected from both types of sur- 
veys. Because the SZE survey has a nearly z-independent 
sensitivity, we find that high-redshift clusters at z > 1 
yield useful constraints, in addition to those studied by 
WS98 in the range < z < 1. Finally, we quantify the 
statistical significance of differences in the models by ap- 
plying a combination of a Kolmogorov-Smirnov (KS) and 
a Poisson test to dN / dzdU,, and obtain constraints using a 
grid of models for a wide range of cosmological parameters. 

This paper is organized as follows. In § we describe 
the main features of the proposed SZE and X-ray surveys 
relevant to this work. In § |^ we briefly summarize our 
modeling methods and assumptions. In § |^, we quantify 
the effect of individual variations of w and of other pa- 
rameters on cluster abundance and evolution. In § |, we 
obtain the constraints on these parameters by considering 
a grid of different cosmological models. In § ra, we discuss 
our results and the implications of this work. Finally, in 
§ 1^, we summarize our conclusions. 

2. CLUSTER SURVEYS 

The observational samples available for studies of clus- 
ter abundance evolution will improve enormously over the 
coming decade. The pres ent samples of tens of intermedi 



ate re dshift clusters (e.g., Gioia et al. 199C , Vikhlinin et al 



1998) will be replaced by samples of thousands of interme- 
diate redshift and hundreds of high redshift (z > 1) clus- 
ters. At a minimum, the analysis of the European Space 
Agency X-ray Multi-mirror Mission (XMM) archive for 
serendipitously detected clusters will yield hundreds, and 
perhaps thousands of new clusters with emission weighte d 
mean temperature measurements ( Romcr ct al. 20001 ). 
Dedicated X-ray and SZE surveys could likely surpass the 
XMM sample in areal coverage, number of detected clus- 
ters or redshift depth. The imminent improvement of dis- 
tant cluster data motivates us to estimate the cosmological 
power of these future surveys. Note that in practice, the 
only survey details we utilize in our analyses are the virial 
mass of the least massive, detectable cluster (as a function 



of redshift and cosmological parameters), and the solid an- 
gle of the survey. We include here a brief description of 
two representative surveys. 

2.1. A Sunyaev-Zel'dovich Effect Survey 

The SZE survey we con sider is that proposed b y Carl- 
strom and collaborators ( Carlstrom et al. 1999 ). This 
interferometric survey is particularly promising, because 
it will detect clusters more massive than ~ 2 x IQ^'^Mq, 
nearly independent of their redshift. Combined, this low 
mass threshold and its redshift independence produce a 
cluster sample which extends, depending on cosmology, to 
redshifts z ~ 3. The proposed survey will cover 12 deg^ 
in a year; it will be carried out using ten 2.5 m telescopes 
and an 8 GHz b andwidth digital c orrelator operating at 
cm wavelengths (Mohr et al. 199!;). The detection limit 
as a function of redshift and cosmology Mi„in(z, Vtrm h) for 
this survey has been studi ed using mock obser vations of 
simulated galaxy clusters ( Holder et al. 2000 ), and we 
draw on those results here. 

Optical and near infrared foUowup observations will be 
required to determine the redshifts of SZE clusters. Given 
the relatively small solid angle of the survey, it will be 
straightforward to obtain deep, multiband imaging. We 
expect that the spectroscopic foUowup will require access 
to a multiobject spectrograph on a 10 m class telescope. 
The ongoing development of infrared spectrographs may 
greatly enhance our ability to effectively measure redshifts 
for the most distant clusters detected in the SZE survey. 

2.2. A Deep, Large Solid Angle X-ray Survey 

We also consider the cosmological sensitivity of a large 
solid angle, deep X-ray imaging survey. The character- 
istics of our survey are similar to those of a proposed 
Small Explorer class mission, called the Cosmology Ex- 
plorer, spear-headed by G. Ricker and D. Lamb. The 
survey depth is 3.6 x 10^ cm^s at 1.5 keV, and the cov- 
erage is IC* deg^ (approximately half the available un- 
obscured sky). We assume that the imaging character- 
istics of the survey are sufficient to allow separation of 
the 10% clusters from the 90% AGNs and galactic stars. 
We focus on clusters which produce 500 detected source 
counts in the 0.5:6.0 keV band, sufficient to reliably esti- 
mate the emission weighted mean temperature in a survey 
of this depth (the external and internal backgrounds sum 
to ~ 1.4 cts/arcmin^). 

To compute the number of photons detected from a 
cluster of a particular flu x, we assume the clusters emit 
Raymond-Smith spectra (Raymond & Smith 1977) with 



i solar abundance, and we model the effects of Galac- 
tic absorption using a constant column density of nn = 
4 X 10^" cm~^. The metallicity and Galactic absorption 
we've chosen are representative for a cluster studied in re- 
gions of high Galactic latitude; when analyzing a real clus- 
ter one would, of course, use the Galactic nn appropriate 
at the location of the cluster. Cluster metallicities vary, 
but for the 0.5:6 keV band, line emission contributes very 
little flux for clusters with temperatures above 2 keV. For 
example, if the cluster metallicity were doubled to | so- 
lar, the conversion between flux and the observed counts in 
the 0.5:6 keV band for this particular survey would vary by 
~ 1.4% and ~ 0.1% for Raymond-Smith spectral models 
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with temperatures kT = 2 keV and 10 keV, respectively. 
We assume that the detectors have a quantum efhciency 



simil ar to the A CIS detectors ( Bautz et al . 1998 
tas et 



Char- 



al. 1998) on the Chandra X-ray Observatory, and 



the energy dependence of the mirror effecti ve area mimics 



that of the mirror modules on ABRIXAS (Friedrich et al 
19980 

The X-ray survey could be combined with the Sloan Dig- 
ital Sky Survey (SDSS) to obtain redshifts for the clusters 
- the redshift distribution of the clusters which produce 
500 photons in the survey described above is well sampled 
at the SDSS photometric redshift limit. 

2.3. Determining the Survey Limiting Mass Afmin 

For our analysis, the most important aspect of both sur- 
veys is the limiting halo mass Mminiz, Om, w^h), as a func- 
tion of redshift and cosmological parameters. More specif- 
ically, we seek the relation between the detection limit of 
the survey a,Tid the correspondinty limitinp; "viria.l mass". 



In oui l modeling below, we will be using the mass function 
of dar k halos obtained in la rge scale cosmological simula- 
tions ( Jenkins et al. 2000| ). In these simulations, halos 
are identified as those regions whose mean spherical over- 
density exceeds the fixed value Sp/pb = 180 (with respect 
to the background density pb, and irrespective of cosmol- 
ogy; see discussion below). In what follows, we adopt the 
same definition for the mass of dark halos associated with 
galaxy clusters. 

In the X-ray survey, Mj^in follows from the cluster X- 
ray luminosity - virial mass relation and the details of the 
survey. We adopt the relation between virial mass and 
temperature obtained in hydrodynamical simulations by 
Bryan & Norman (1998), 



7-3/2 



(1) 



where H{z) = HqE{z) is the Hubble parameter at redshift 
z, a = 1.08 is a normalization determined from the hydro- 
dynamical simulations, and Ac is the enclosed overdensity 
(relative to the critical density) which defines the cluster 
virial region. The normalization a is found to be relatively 
insensitive to cosmological parameters, and the redshift 
evolution of Equation |lj appears to be consistent with the 
hydrodynami cal simulations in those models where it has 
been tested ( Bryan fc Norman 1998 ). Here we assume 
that Equation"! holds in all cosmologies with the same 
value of a (see §6.2 for a discussion of the effects of er- 
rors in the mass-temperature relation) , and use the fitting 
formulae for Ac provided by WS98, which includes the 
case w ^ —1. Finally, we convert Mvir from Equation |] to 
the mass Migo enclosed within the spherical overdensity of 
6p/p = 180 (with respect to the background density), as- 
suming that the halo profile is well described by the NFW 
model with concentration c — 5 (Navarro, Frenk & White 
1997, hereafter NFW). 

We next utilize Equation ||, together with the relation 
between bolometric luminosity and temperature found by 
Arnaud & Evrard (1999), to find the limiting mass of a 
cluster that produces 500 photons in the 0.5:6.0 keV band 
in a survey exposure. For these calculations we assume 
that the luminosity-temperature relation does not evolve 



with red shift, consistent with the cu rrently available obser- 
vations ( Mushotzky fc Scharf 1997 ; relaxing this assump- 
tion is discussed below in § 

For an interferometric SZE survey, the relevant observ- 
able is the cluster visibility V, which is the Fourier trans- 
form of the cluster SZE brightness distribution on the sky 
as seen by the interferometer. The visibility is proportional 
to the total SZE flux decrement S'y, 



V oc S^{M, z) (X ficM 



d\{z) 



(2) 



where (T'e)„ is the electron density weighted mean tem- 
perature, M is the virial mass, ficM is the intracluster 
medium mass fraction and is the angular diameter dis- 
tance. We normalize this relation usi ng mock observations 
of numerical cluster simulations (sec Mohr & Evrard 1997 



and Mohr, Mathiesen & Evrard 1999) carried out in three 
different cosmological models, including noise ch aracteris- 
tics appr opriate to the proposed SZE array (see Holder et 



al. 200C for more details). The ICM mass fraction is set to 
fiCM = 0.12 in all three cosmological models. This mass 
fraction is consistent with analyses of X-ray emission from 
well defined samples if i?o — 65 km s^^ Mpc^^, our fidu- 
cial value. Note that we use the same ficM = 0.12 in all 
our cosmological models rather than varying it with the 
iJo scaling appropriate for analyses of cluster X-ray emis- 
sion. In the discussion which follows, this choice allows us 
to focus solely on the cosmological discriminatory power 
of cluster surveys; naturally, in interpreting a real cluster 
survey one would likely allow ficAi to vary with Hq. 




Fig. 1. — Limiting cluster virial masses (Miso) for detection in the 
X-ray survey (upper pair of curves) and in the SZE survey (lower 
pair of curves). The solid curves show the mass limit in our fiducial 
fiat ACDM model, with w = -1, = 0.3, and h = 0.65, and 
the dotted curves show the masses in the same model except with 
w = -0.5. 

Note that for a fiux limited survey, the limiting mass 
in equation || is sensitive to cosmology through its depen- 
dence on dA and the definition of the virial mass M. We 
adopt the simulation-normalized value of M^-^^-^{z) in our 
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fiducial cosmology as a template, and then we rescale this 
relation to determine M„iin(z) in the model of interest us- 
ing the relation 



h* 



h*d\{z) 



6/5 



(3) 



Here the superscript * refers to quantities in the ACDM 
reference cosmology, and we have used the scaling of virial 
mass with temperature (Eqn. |l|): M cx {T^fJ^. We tested 
this scaling by comparing it to mock observations in simu- 
lations of two different cosmologies (open CDM and stan- 
dard CDM), and found that agreement was better than 
~ 10% in the redshift range < z < 3. Finally, in the nu- 
merical simulations used to calibrate Equation the halo 
mass was defined to be the total mass enclosed within a 
region whose mean spherical interior density is 200 times 
the critical density. As in the X~ray case, we convert 
Afmin(^) from Equation || to the desired mass Migo by 
assuming that the halo profile follows NEW with concen- 
tration c = 5. 

The mass limits we derive for both surveys are shown in 
the redshift rang e < z < 3 in Figure both for ACDM 
and for aw — —0.5 universe. The SZE mass limit is nearly 
independent of redshift, and changes little with cosmology. 
As a result, the cluster sample can extend to z « 3. In 
comparison, the X-ray mass limit is a stronger function of 
w, and it rises rapidly with redshift. For the X-ray survey 
considered here the number of detected clusters beyond 
z « 1 is negligible. 

These mass limits incorporate some simplifying assump- 
tions that have not been tested in detail (although we con- 
sider small variations of the mass limits below) . Our goal is 
to capture the scaling with cosmological parameters and 
redshift as best as presently possible. However, we em- 
phasize that further theoretical studies of the sensitivities 
of these scalings to, for example, energy injection during 
galaxy formation will be critical to interpreting the survey 
data. In the case of the X-ray survey, the cluster sam- 
ple will have measured temperatures, allowing the limiting 
mass to be estimated independent of the cluster luminos- 
ity. In the case of the SZE survey, deep X-ray foUowup 
or multifrequency SZE foUowup observations should yield 
direct measurements of the limiting mass. 

3. ESTIMATING THE CLUSTER SURVEY YIELD 

To derive cosmological constraints from the observed 
number and redshift distribution of galaxy clusters, the 
fundamental quantity we need to predict is the comov- 
i ng cluster mass function . The Press-Schechter formalism 
( [Press fc Schechter 1974 ; hereafter PS), which directly pre- 
dicts this quantity in any cosmology, has been shown to 



be in i rcaaonably good agreement (i.e. — to within a factor 



of two ) with roeulte of N-body eimulatione, in coe mologioe 
and halo mass ranges where it has been tested ( [ Lacey fc 



Cole 1994, Gross et al 1998, Lee k Shandarin 1999 



merical simulations have only recently reached the large 
size required to accurately determine the mass function of 
the rarest, most massive objects, such as galaxy clusters 
with M > 10^^ Mq. 

In this paper, we adopt the halo mass function found in 
a series of recent large-scale cosmological simulations by 
Jenkins et al. 2000. The results of these simulations are 



particularly well-suited for the present application. The 
large simulated volumes allow a statistically accurate de- 
termination of the halo mass function; for halo masses of 
interest here, to better than <^ 30%. In addition, the 
mass function is computed in three different cosmologies 
at a range of redshifts, and found to obey a simple 'univer- 
sal' fitting formula. Although this does not guarantee that 
the same scaling holds in other, untested cosmologies, we 
make this simplifying assumption in the present paper. In 
the future, the validity of this assumption has to be tested 
by studying the numerical mass function across a wider 
range of cosmologies. 

Generally, the simulation mass function predicts a sig- 
nificantly larger abundance of massive clusters than does 
the PS formula. For sake of definiteness, we note that 
in the simulations, halos are identified as those regions 
whose mean spherical overdensity exceeds the fixed value 
Sp/pb = 180 with respect to the background density pi,. 
This is somewhat different from the typical halo definition 
within the context of the PS formalism, where the over- 
density, relative to the critical density, is taken to be that 
of a collapsing spherical top-hat at virialization. 

Following Jenkins et al. 2000, we assume that the co- 
moving number density [dn / dM)dM of clusters at redshift 
z with mass M ± dM /2 is given by the formula, 



— -(z, M) — 0.315 — — - 

dM^ ' M au dM 



exp 



- 10.61 - log{D,aM. 



|3.8 



(4) 

where aM is the r.m.s. density fluctuation, computed on 
mass-scale Al from the present-day linear power spectrum 
( Eisenstein fc Hu 1998 ), Dz is the linear growth function, 
and po is the present-day mass density. The directly ob- 
servable quantity, i.e. the average number of clusters with 
mass above Mmin at redshift z ± dz/2 observed in a solid 
angle dQ is then simply given by 



dN 
dzdfl 



(^) 



dV 
dzdfl 



(^) 



dM 



dn 
dM 



(5) 



where dV/dzdQ is the cosmological volume element, and 
Mminiz) is the limiting mass as discussed in section 2.3 



Equations ^ and |^ depend on the cosmological param 
eters through po, D^, and dV/dzdil, in addition to the 
mild dependence of ctm on these parameters through the 
power spectrum (although the dependence on the power- 
spectrum is more pronounced in the X-ray survey, where 
the limiting mass varies strongly with redshift). Note that 
the comoving abundance dn/dM is exponentially sensitive 
to the growth function D^. We use convenient expres- 
sions for dV/dzdVl and i n open and flat Vl a cosmologies 
available in the literature (iPeebles 198(]| , [CarroU, Press &| 
Turner 1992| , [Eisenstein 1996| ). In the case of cosmologies 



with K7 7^ — 1, we have evaluated dV/dzdfl numerically, 
but used the fitting formulae for obtained by WS98, 
which are accurate to better than 0.3% for the cases of 
constant w's considered here. 

3.1. Normalizing to Local Cluster Abundance 

To compute dN/ dzdfl from equation we must choose 
a normalization for the density fluctuations ctm- This is 
commonly expressed by erg; the present epoch, linearly 
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extrapolated rms variation in the density field filtered 
on scales of 8h~^ Mpc. To be consistent in our analy- 
sis, we choose the normalization for each cosmology by 
fixing the local cluster abundance above a given mass 
Mnm = 10^'^h~^ Mq. In all models considered, we set the 
local abundance to be 1.03 x 10"^ {h/0.65)^ Mpc'^, the 
value derived in our fiducial ACDM model (see below). We 
have chosen to normalize using the local cluster abundance 
(upto a factor h^) above mass M^m rather than above a 
particular emission weighted mean temperature fcTmn, be- 
cause this removes the somewhat uncertain cosmological 
sensitivity of the virial mass temperature {M—T^) relation 
from the normalization process; spherical tophat calcula- 
tions suggest a significant offset in the M — normaliza- 
tion of the open and flat = 0.3 models whic h hydro 
dynar nical simulations do not seem to reproduce ( Evrard 
Metzl hr fc Navarro 1996| , [Bryan fc Norman 199^ , [Viana 
Liddl( t 19991 )^ 

An alternative approach to the above is to regard dg 
a "free-parameter", on equal footing with the other pa- 
rameters we let float below. This possibility will be dis- 
cussed further in § 6. Here we note that our normaliza- 
tion approach is sensible, because the number density of 
nearby clusters can be measured to within a factor of , 
and the masses of nearby clusters can be measured di- 
rectly through several independent means; these include 
the assumption of hydrostatic equilibrium and using X-ray 
images and intracluster medium (ICM) temperature pro- 
flies, weak lensing, or galaxy dynamical mass estimates. 
The only cosmological sensitivity of these mass estima- 
tors is their dependence on the Hubble parameter h] we 
include this h dependence when normalizing our cosmo- 
logical models. Note that previous derivations of tJs (e.g. 
Viana & Liddle 1993; Pen 1998) in various cosmologies 
from the local cluster abundance N{> kT) above a flxed 
threshold temperature fcTmin 7keV yielded a constraint 



with the approximate scaling crgf^ 



1/2 



0.5. We find a 



similar relation when varying away from our fidu- 
cial cosmology; however, we note that if a ~ 5 times 
smaller threshold temperature were used, the constrained 
combination would be quite different, trs^^m ^constant. 
Since our adopted normalization is based on mass, rather 
than temperature, in general, we find still different scal- 
ings. As an example, when h = 0.65 and w = —1 are 
kept fixed, our normalization procedure translates into 
a8(f]„/0.3)0-85 « 0.9. 

3.2. Fiducial Cosmological Model 

The parameters we choose for of our fiducial cosmolog- 
ical model are (SIa, fim, /i, cg, n) — (0.7,0.3,0.65,0.9,1). 
This fiat ACDM model is motivated as a "best-fit" model 
that produces a local cluster abun dance consistent with 
observations (Viana & Liddle 1999), an d satisfies the cur- 
rent con straints from CMB anisotropy dLange et al. 2000| , 
see also White, Scott & Pierpaoli 2000), high-z SNe, and 
other observations (Bahcall et al. 1999). We have as- 
sumed a baryon density of Q.i,h? — 0.02, consistent with 
recent D/H measurements (e.g. Buries & Tytler 1998). 
Note that the power spectrum index n is not important 
for the analysis presented here, because we normalize on 
cluster scales erg , and we find that this minimizes the effect 
of varying n on the density fluctuations relevant to cluster 



formation. 



4. EXPLORING THE COSMOLOGICAL SENSITIVITY 

In this section, we describe how variations of the indi- 
vidual parameters w, and h, as well as the cosmological 
dependence of the limiting mass Mmin, affect the cluster 
abundance and redshift distribution. This will be useful 
in understanding the results of the next section, when a 
full grid of different cosmologies is considered. We then 
describe our method of quantifying the statistical signifl- 
cance of differences between the distributions dN/dz in a 
pair of different cosmologies. 



4.1. Single Parameter Variations 

The surface density of clusters more massive than Mmin 
depends on the assumed cosmology mainly through the 
growth function D{z) and volume element dV/dzdfl, as 
well as through the cosmology-dependence of the limit- 
ing mass Mmin itself. In the approach described in sec- 
tion 1^, once a cosmology is specified, the normalization 
of the power spectrum erg is found by keeping the abun- 
dance of clusters at z = constant. We therefore con- 
sider only three "free" parameters, w, h, 17^, specifying 
the cosmology. We assume the universe to be either fiat 
{Qq = 1 — Qm), or open with Qq = 0. 
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Fig. 2. — Effect of changing Qm when all other parameters are 
held fixed. The four panels show (clockwise from upper left) the 
surface density of clusters at redshift z; the linear growth function; 
the volume element in units of Mpc^ sr~^ redshift"^; and the co- 
moving cluster abundance. The solid curve shows our fiducial flat 
ACDM model, with w = —1, Qrn = 0.3, and h = 0.65. Also shown 
are models with Q = 0.27 (dotted curve); Q = 0.33 (short-dashed 
curve); and OCDM models with ft = 0.27,0.30,0.33 (long-dashed 
curves, top to bottom). 
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Fig. 3. — Effect of changing w when all other parameters are held 
fixed. The solid curve shows our fiducial flat ACDM model, with 
w = —I, Qrn = 0.3, and h = 0.65. The dotted curve is the same 
model with to = —0.6, the short-dashed curve with w = —0.2, and 
the long-dashod curve is an open CDM model with Qrn = 0.3. 

4.1.1. Changing flm 

The effects of changing $7,,,, are demonstrated in Fig- 
ure ||. The curves correspond to a flat ACDM universe 
with {h = 0.65, w = -1), and n„i = 0.27 (dotted), 
— 0.30 (sohd), and flm — 0.33 (short-dashed). In ad- 
dition, the long-dashed curves show the same three mod- 
els (top to bottom), assuming open CDM with JIa = 0. 
The top left panel shows the total number of clusters in 
a 12 square degree field, detectable down to the constant 
SZE decrement S'min- As discussed in section 2.3 above, a 
constant S'min implies a redshift and cosmology-dependent 
limiting mass Mmin- In the SZE case, we find that if we 
had not included this effect, the sensitivity to would 
have been somewhat stronger. Several conclusions can be 
drawn from Figure||. Overall, the top left panel shows that 
a decrease in J7„j increases the number of clusters (and vice 
versa) at all redshifts. Note that the dependence is strong, 
for instance, a 10% decrease in increases the total num- 
ber of clusters by 30% in either ACDM or OCDM cos- 
mologies. As emphasized by Bahcall & Fan (1998), Viana 
& Liddle (1999) and others, this makes it possible to es- 
timate an upper limit on using current, sparse data 
on cluster abundances (i.e. only a few high-z clusters). 
A second important feature seen in the top left panel is 
that the shape of the redshift distribution is not changed 
significantly, a conclusion that holds both in ACDM and 
OCDM. Finally, the remaining three panels reveal that the 
effects of rim arise mainly from the changes in the comov- 
ing abundance (bottom left panel). In flat ACDM, flm has 
relatively little effect on the volume or the growth func- 
tion, and the comoving abundance is determined by the 
value of (Tg that keeps the the local abundance constant at 
z — (we find tT8=0-83 for fi^ = 0.33 and as = 1-00 for 
rim = 0.27). In addition, we find that the change in the 
shape of the underlying power spectrum with flm enhances 



the differences caused by Vim (when we artificially keep the 
power spectrum at its flm — 0.3 shape, we find cr8=0.84 
for rim = 0.33). We also note that the volume element and 
the comoving abundance act in the same direction: a lower 
Qm increases both the comoving abundance and the vol- 
ume element. In OCDM, the growth function has a larger 
effect, and relative to ACDM, the redshift distribution is 
much flatter. 
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Fig. 4. — Effect of changing h when all other parameters are held 
fixed. The ACDM model of Figure H is shown (solid curve) together 
with models with h = 0.55 (dotted curve); h = 0.80 (short-dashed 
curve); and OCDM models with h = 0.55,0.65,0.80 (long-dashed 
curves, top to bottom). 

4.1.2. Changing w 

The effects of changing w are demonstrated in Figure |[ 
The figure shows models with {Um — 0.3, h — 0.65) and 
with three different w's: w — —1 (solid curve), w — —0.6 
(dotted curve), and w = —0.2 (short-dashed curve). In 
addition, we show the result from an open CDM model 
with (il = 0.3, /i = 0.65; long-dashed curve). The fig- 
ure reveals that increasing w above w = —1 causes the 
slope of the redshift distribution above z w 0.5 to flatten, 
increasing the number of high-z clusters. Furthermore, 
"opening" the universe has an effect similar to increas- 
ing w. The other three panels demonstrate the reason for 
these scalings. The top right panel shows that the growth 
function is flatter in higher w models, signiflcantly increas- 
ing the comoving number density of high-redshift clusters 
(bottom left panel). The volume element (bottom right 
panel) has the opposite behavior, in the sense the volume 
in higher-w models is smaller, which tends to balance the 
increase in the comoving abundance caused by the growth 
function in the range < z <i 0.5; but for higher redshifts, 
the growth function "wins" . An important conclusion seen 
from Figure || is that both the total number of clusters 
as well as the shape of their redshift distribution, signifi- 
cantly depends on w. We also note that in the SZE case, 
our sensitivity to w has been enhanced by the cosmological 
dependence of the mass limit (opposite to what we found 
for the Om-sensitivity, which we found was weakened by 
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the same effect). 

4.1.3. Changing h 

Figure ^ demonstrates the effects of changmg h. Three 
ACDM models are shown with (f7„i = 0.30, = —1), 
and h = 0.55 (dotted curve), h — 0.65 (sohd curve), and 
h = 0.80 (short -dashed curves). The long-dashed curves 
correspond to OCDM models with the same parameters 
(top to bottom). Comparing the top right panel with 
that of Figure g, the qualitative behavior of dN/dz under 
changes in h and flm are similar: decreasing h increases 
the total number of clusters, but does not considerably 
change their redshift distribution. However, the sensitiv- 
ity to h is significantly less: the total number of clusters 
is seen to increase by ~ 25% only when h is decreased 
by the same percentage. Note that the growth function is 
not effected by h, and the h sensitivity is driven by our 
norm aliza tion process, which fixes the abundance at z = 
(see § 3.1). Since the volume scales as oc h~^, we fix the co- 
moving abundance to be proportional to oc /i^. As a result, 
dN/dzdil is nearly independent of h. In fact, the entire h- 
dependence is attributable to the small change caused by 
h in the shape of the power spectrum (for a pure power- 
law spectrum, there would be no /i-dependence, and the 
three curves for the flat universe in the top left panel of 
Figure ^ would look identical) . 
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Fig. 5. — Effect of changing ui (upper panels) or Qrn (lower pan- 
els) when all other parameters are held fixed, including the mass 
limit. The types of the curves corresnondJjO the different models in 
the SZE survey, as shown in Figures Q & y. 

4.1.4. Abundances in the X-ray Survey 

The evolution of the cluster abundance, and its sensi- 
tivity to ri„j and w in the X-ray survey are shown in Fig- 
ure |6|. Because of the much larger solid angle surveyed, the 
numbers of clusters is significantly larger than in the SZE 
case, despite the higher limiting mass (cf. Fig |l|). Nev- 
ertheless, the general trends that can be identified in the 
X-ray sample are similar to those in the SZE case. Rais- 
ing w increases the total number of clusters, and flattens 



their redshift distribution. As in the SZE survey, raising 
rim decreases the total number of clusters. 

4.2. Effects of the Limiting Mass Function 

Finally, we examine the extent to which the above 
conclusions depend on the cosmology and redshift- 
dependence of the limiting mass M^i^. 
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z 

Fig. 6. — Effect of changing w (upper panels) or f2m (lower pan- 
els) when all other parameters are held fixed in the X-ray survey. 
Note the much larger numbers of clusters in comparison to the SZE 
survey. In the top panel, the curves correspond to to = —1 (solid), 
UI = —0.6 (dotted) and w = —0.2 (dashed). In the bottom panel, 
the curves correspond to = 0.3 (solid), = 0.27 (dotted) and 
f^m = 0.33 (dashed). 



4.2.1. The SZE Survey 

We first compute cluster abundances above the fixed 
mass Mmin = lO^^/i~^M0, characteristic of the SZE sur- 
vey detection threshold in the range of cosmologies and 
rcdshifts considered here. The results are shown in Fig- 
ure ^: the bottom panels show the surface density and 
comoving abundance when f2m is changed (the models 
are the same as in Figure ||), and the top panels show 
the same quantities under changes in w (the cosmological 
models are the same as in Figure A comparison be- 
tween Figures ^ and ^ gives an idea of the importance of 
the mass limit. The general trend seen in Figure ^ remains 
true, i.e. increasing w flattens the redshift distribution at 
high-z. However, when a constant M^^^^ is assumed, the 
"pivot point" moves to slightly higher redshift, and the to- 
tal number of clusters becomes less sensitive to w. Similar 
conclusions can be drawn from a comparison of Figure |^ 
with the bottom two panels of Figure 0: under changes in 
fim the general trends are once again similar, but the dif- 
ferences between the different models are amplified when 
a constant Mmin is used. In summary, we conclude that 
in the SZE case (1) the variation of the mass limit with 
redshift and cosmology has a secondary importance, and 
(2) it weakens the r2„i dependence, but strengthens the w 
dependence. 
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4.2.2. The X-ray Survey 

In comparison to the SZE survey, the X-ray mass hmit 
is not only higher, but is also significantly more dependent 
on cosmology (cf . Fig |l|) . On the other hand, the X-ray 
sample goes out only to the relatively low redshift z = 1, 
where the growth functions in the different cosmologies di- 
verge relatively little. This suggests that in the X-ray case 
the mass limit is more important than in the SZE survey. 
In order to separate the effects of the changing mass limit 
from the change in the growth function and the volume 
element, in Figure |^ we show the sensitivity of dN/dz to 
changes in flm and w, without including the effects from 
the mass limit. The same models are shown as in Figure ^ 
except we have artificially kept the mass limit at its value 
in the fiducial cosmology. The figure reveals that essen- 
tially all of the w-sensitivity seen in Figure ^ is caused 
by the changing mass limit; when M^nin is kept fixed, the 
cluster abundances change very little. On the other hand, 
comparing the bottom panels of Figures ^ and |^ shows that 
including the scaling of the mass limit somewhat reduces 
the rim dependence, just as in the SZE case. 




s 



-9 



0.3 0.4 0.6 0.8 1 
Z 



O.S 0.4 0.6 o.e 

z 



Fig. 7. — Effect of changing w (upper panels) or Qm (lower panels) 
when all other parameters are held fixed in an X— ray survey, and 
the survey mass limit is held fixed at its fiducial value, irrespective 
of cosmology. A comparison with Figure ra shows that nearly all of 
the ui-sensitivity is accounted for by the cosmology— dependence of 
the limiting mass. On the other hand, the Qm-sensitivity is caused 
mostly by the growth function. 



4.3. Overview of Cosmological Sensitivity 

In summary, we conclude that changes in w modify both 
the normalization and the shape of the redshift distribu- 
tion of clusters, while changes in n„i or h effect essentially 
only the overall amplitude. This suggests that changes in 
w can not be fully degenerate with changes in either 0^ 
or h (or a combination), making it possible to measure w 
from cluster abundances alone. These conclusions hold ei- 
ther for clusters above a fixed detection threshold in and 
SZE or X-ray survey, or for a sample of clusters above 



a fixed mass. We find that the sensitivity to flm arises 
mostly through the growth function, both in the SZE and 
X-ray surveys. This sensitivity is slightly weakened by the 
scaling of the limiting mass M^i^ with flm- We find that 
the w sensitivity is also dominated by the growth func- 
tion in the SZE survey, which goes out to relatively high 
redshifts; but the sensitivity to w is enhanced by the w- 
dependence of M^i^- In comparison, in the X-ray survey, 
which only probes relatively low redshifts, nearly all of 
the w-sensitivity is caused by the cosmology-dependence 
of the limiting mass, rather than the growth function. 

5. CONSTRAINTS ON COSMOLOGICAL PARAMETERS 

We derive cosmological constraints by considering a 3- 
dimensional grid of models in fi™, h, and w. As described 
above, we first find erg in each model, so that all models are 
normalized to produce the same local cluster abundance 
at 2 = 0. We then compute dN/dzdO. in these models for 
0.2 < f7™ < 0.5, 0.5 < /i < 0.9, and -1 < w < -0.2. 
The range for w correspon ds to that allowed by current 
astrophysical observations (Wang et al. 200C); although 
recent observations of Type la SNe suggest the stronge r 
constraint w <, — 0.6 ( Pcrlmutter, Turner fc White 1999 ). 



5.1. Comparing dN/dz in Two Different Cosmologies 

The main goal of this paper is to quantify the accuracy 
to which w can be measured in future SZE and X-ray 
surveys. To do this, we must answer the following ques- 
tion: given a hypothetical sample of iVtot clusters (with 
measured redshifts) obeying the distribution dNA/dz of 
the test model (A) cosmology, what is the probability 
Ptot(^, -S) that the same sample of clusters is detected 
in the fiducial (B) cosmo logy , with distribution dNs/dzl 
We have seen in section 4.1 that the overall amplitude, 
and the shape of dN/dz are both important. Motivated 
by this, we define 



Ptot{A,B) = Po{A,B) X PM,B) 



(6) 



where P{){A, B) is the probability of detecting iV^i^tot clus- 
ters when the mean number is NB,tot, and Pz^A, B) is the 
probability of measuring the redshift distribution of model 
(A) if the true parent distribution is that of model (B). We 
assume Pq is given by the Poisson distribution, and we use 
the Kolmogorov-S mirnov (KS) test to compute Pz{A,B) 
(Press et al. 1992). The main advantage of this approach, 
when compared to the usual tests, is that we do not 
need to bin the data in redshift. 

For reference, it is useful to quote here some examples 
for the probabilities, taking (ri„i = 0.3, h — 0.65, w = — 1) 
as the fiducial (B) model. For example, closest to this 
model in Figure H is the one with w — —0.6. For this 
case, we find Pq — 0.25 and Pz — 0.1 for a total proba- 
bility of Ptot — 0.025. In other words, the two cosmolo- 
gies could be distinguished at a likehhood of 1.2cr using 
only the total number of clusters, at 1.6(T using only the 
shape of the redshift distribution, and at the 2.3(7 level 
using both pieces of information. In this case, the distinc- 
tion is made primarily by the different redshift distribu- 
tions, rather than the total number of detected clusters. 
Taking the = 0.33 ACDM cosmology from Figure | 
as another example for model (A), we find Po — 0.0075 
(=2.7ct), Pz = 0.78 (=0.3cr), and a total probability of 
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Fig. 8. — Contours of 1, 2, and 3a likelihood for different models when they are compared to a fiducial flat ACDM model with Qrn = 0.3 
and h = 0.65, using the SZE survey. The three panels show three different cross— sections of constant total probability at fixed values of h 
(0.55,0.65, and 0.80) in the investigated 3-dimensional Q,m,w,h parameter space. 

^tot = 0.0058 (=2.8(t). Not surprisingly, the shape of the 
redshift distribution does not add significantly to the sta- 
tistical difference between these two models, which differ 
primarily by the total number of clusters. 

5.2. Expectations from the Sunyaev-Zel'dovich Survey 

Figure || shows contours of f , 2, and 3a for the total 
probability Ptot for models when compared to the fiducial 
flat ACDM model. For reference, we note that the total 
number of clusters in the SZE survey in our fiducial model 
is Ri 100, located between < z < 3. The three panels 
show three different cross-sections of the investigated 3- 
dimensional U,m,,h,w parameter space, taken at constant 
values oi h = 0.55, 0.65, and 0.80, spanning the range of 
values preferred by other observations. The most striking 
feature in this figure is the direction of the contours, which 
turn upwards in the w, Vim plane, and become narrower for 
larger values of w. We find that the trough of maximum 
probability for fixed h = 0.65 is well described by 



(r!„, -0.3)(u; + l)-5/2 = 0.1, 



(7) 



with further constant shifts in caused by changing h. 
The ±3(7 width enclosed by the contours around this re- 
lation is relatively narrow in (±10%). In a ACDM 
case, even when a large range of values is considered for 
h (0.45 < h < 0.90), the constraint 0.26 <, Vt^ ^ 0.36 
follows; when w ^ —1 is considered, the allowed range 
widens to 0.27 ^ fim ^ 0.41. On the other hand, a wide 
range of w's is seen to be consistent with w = —\: the 
largest value shown, w w —0.2 is approximately 3(t away 
from w = —1, and w = —0.6 is allowed at Icr. Note that 
h is not well determined, i.e. the contours look similar for 
all three values of h, and la models exist for any value 
of h in the range 0.5 ^ h <, 0.9. This is not surprising, 
as Figure ^ shows dN/dzdVl is insensitive to the value of 
/i, with only a mild ^.-dependence through the non-power 
law shape of the power spectrum. 

5.3. Expectations from the X-ray Survey 

The total number of clusters in the X-ray survey in our 
fiducial model is ~ 1000, ten times that in the SZE sur- 
vey; all X-ray clusters are located between < z < 1. 
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Fig. 9. — Contours of 1, 2, and 3cr likelihood for models when they are compared to a fiducial flat ACDM model, as in Figure |8| but for 
the X-ray survey. 
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Figure ^ contains expectations for the X-ray survey; we 
show contours of 1, 2, and Scr probabihties relative to the 
fiducial ACDM model. The qualitative features are similar 
to that in the SZE case, but owing to the larger number of 
clusters, the constraints are significantly stronger and the 
contours are narrower. However, the contours extend fur- 
ther along the w axis, and the largest value of w allowed at 
a probability better than 3a is w > —0.2 (assuming that 
the values of and h are not known) . Although the con- 
tours are narrower than in the SZE case, assuming that h 
and w are unknown, the allowed range of flm is similar to 
that in the SZE case, 0.26 ^ ^ 0.42. Note that be- 
cause of the shape and direction of the likelihood contours, 
a knowledge of h would not significantly improve this con- 
straint (although if h is found to be low, then the lower 
limit in would increase). Finally, assuming that both 
h and flm are known to high accuracy (« 3%), the allowed 
3cr range on w could be reduced to — 1 < it; ^ — 0.85. 

6. RESULTS AND DISCUSSION 

6.1. Total Number vs. the Redshift Distribution 

Our main results are presented in Figures || and ^, which 
show the probabilities of various models relative to a fidu- 
cial ACDM model in the SZE and X-ray surveys. As 
demonstrated by these figures, the cluster data determine 
a combination of f2„i and w. In the absence of external 
constraints on and h, w as large as —0.2 differs from 
w = —1 by 3cr; while w — —0.6 would be la away from our 
fiducial ACDM cosmology. Owing to the larger number of 
clusters in the X-ray survey, the constrained combination 
of fl„i and w is significantly narrower than in the SZE 
survey; the direction of the contours is also somewhat dif- 
ferent. As a result, analysis of the X-ray survey could 
distinguish a w w —0.85 model from ACDM at 3a signifi- 
cance, provided Qm is known to an accuracy of ^ 3% from 
other studies. 

It is interesting to ask whether these constraints arise 
mainly from the total number of detected clusters, or from 
their redshift distribution. To address this issue, in Fig- 
ure ^ we show separate likelihood contours for the prob- 
ability Pq (total number of clusters, left panels), and for 
the probability Pz (shape of redshift distribution, right 
panels). In the SZE case, the contours of likelihood from 
the shape information alone are broad, and adding these 
constraints to the Poisson-probability plays almost no role 
in the range w ^ — 0.7 (the contours of Ptot and Pq are 
very similar). However, at larger w, the shape becomes 
increasingly important. Adding in this information signif- 
icantly reduces the allowed region relative to the Poisson- 
probability alone at w ^ — 0.7. It is the combination of 
the Pq and Pz contours that allows ruling out w ^ — 0.2 
at the 3a level. Note that the difference in shapes arises 
mostly from the high-redshift {z ^ 1) clusters (cf. Fig. 

In the X-ray case (bottom panels in Fig. |l^), the situ- 
ation is different, because the contours of Pq and Pz are 
both much narrower. As a result, the contours for the 
combined likelihood are somewhat reduced, but they still 
reach to w ~ —0.2 (at ^ 2a). Note that as in the SZE sur- 
vey, the redshift distribution (of clusters primarily in the 
< z < 1 range) plays an important role. As Figures ^ 
and ^ show, the total number of clusters can be adjusted 
by changing Vim and h. In terms of the total number of 



clusters, w is therefore degenerate both with and h: 
raising w lowers the total number, but this can always be 
offset by a change in Sim and/or h. The bottom left panel 
in Fig. |lO| reveals that based on Pq alone, w ~ —0.2 (and 
flm = 0.43) can not be distinguished from ACDM even 
at the la level. On the other hand, the middle panel in 
Fig. H shows that when the shape information is added, 
w < — 0.2 follows to 2a significance. 
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Fig. id. — Likelihood contours of 1, 2, and 3(7 probabilites as in 
Figures M and H, but when only the total number of clusters (left 
panels), or only the redshift distributions (right panels) are used to 
compute the likelihoods between two models. 



6.2. Discussion of Possible Systematic Uncertainties 

Our results imply that the cluster abundances in the 
SZE and X-ray surveys can provide useful constraints on 
cosmological parameters, based on statistical differences 
expected among different cosmologies. The purpose of this 
section is to summarize and quantify the various system- 
atic uncertainties that can affect these constraints. 

Knowledge of the Limiting Mass Mmin- Our conclusions 
above are dependent on the chosen limiting mass, which 
is a function of both redshift and cosmology. From the 
discussion in 



4.1 



we have seen that the limiting mass 
plays a secondary role in the SZE survey, where the bulk 
of the constraint comes from the growth function. In com- 
parison, we find that Mmin plays an important role in the 
X-ray survey. To demonstrate the importance of the mass 
limit explicitly, in Figure |ll| we show the likelihood con- 
tours in the f2„i — w plane when the variations of the limit- 
ing mass with cosmology are not taken into account. Not 
surprisingly, this makes the contours somewhat narrower, 
but nearly parallel to w ~ this is consistent with our find- 
ing in Figure ^ that the mass limit accounts for nearly all 
of the w-dependence, but it reduces the flm dependence. 
Figure ^ demonstrates the need to accurately know the 
limiting mass Mmin, and its cosmological scaling, in the 
X-ray survey. 
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Fig. 11. — Likelihood contours for a RxecLh = 0.65 in the X- 
ray survey, as in the middle panel of Figure (A but zooming in for 
clarity. The added (nearly horizontal) contours shows the allowed 
region when variations of the limiting mass with cosmology are not 
taken into account. 

Because our proposed cluster sample will have nieasured 
X-ray temperatures, the uncertainty in our knowledge of 
the limiting mass will likely be dominated by the theo- 
retical uncertainties of the M — T relation. In order to 
quantify the effect of such errors, we have performed a set 
of simple modifications to our modeling of the constraints 
from the X~ray survey. In all cases, we adopt the same 
M — T relations as we did before (cf. eq. |l]). However, in 
the fiducial model, we use a limiting mass that is altered 
by either ±5% or ±10% from the mass inferred from this 
M — T relation. This mimics a situation where the theo- 
retical M — T relation we apply is either 5% or 10% away 
from the relation in the real universe. In a second set of 
calculations, we mimic a situation where the slope of the 
M — T relation is incorrectly modeled; i.e. we alter this 
slope in the fiducial model to a = 1.5 ± 0.05. The devi- 
ations to the likelihood contours caused by these offsets 
are demonstrated in Figure |l2|, which shows the effects of 
the offset in the M — T normalization, and in Figure O, 
which shows the effects of the offsets in the slope. As 
the figures reveal, the contours shift relatively little under 
these changes. We conclude that the results we derive are 
robust, as long as we can predict the M — T relation to 
within - 10%. 

In our approach, we have attempted to utilize the whole 
observed cluster sample, down to the detection threshold: 
we had to therefore include the above cosmological depen- 
dencies. In principle, measured cluster velocity dispersions 
and X-ray temperatures (both of which are cosmology in- 
dependent) could be utilized to improve the constraints, 
i.e. by selecting sub-samples that maximize the differ- 
ences between models. Further work is needed to clarify 
the feasibility of this approach, as well as to quantify the 
accuracy to which the dependence of Mmin on Q.m, h, w, 
and z can be predicted. 

Evolution of Internal Cluster Structure. Further work is 



also required to test the cluster structural evolution models 
we use. For the X-ray survey, we have assumed that the 
cluster luminosity-temperature relat ion does not evolve, 
consi stent with current observations ( Mushotzky & Scharl 
1997D , and in the SZE survey, we have adopted the struc- 



tural evolution found in state of the art hydrodynamical 
simulations. Because of the sensitivity of the survey yields 
to the limiting mass, cluster structural evolution which 
changes the observability of high redshift clusters can in- 
troduce systematic errors in cosmological constraints: for 
example, both low flm cosmologies and positive evolution 
of the cluster luminosity-temperature relation increase the 
cluster yield in an X-ray survey. SZE surveys are generally 
less sensitive to evolution than X-ray surveys, because the 
X-ray luminosity is heavily dependent on the core struc- 
ture (e.g., the presence or absence of cooling instabilities), 
whereas the SZE visibility depends on the integral of the 
ICM pressure over the entire cluster (Eqn. ||). We are test- 
ing these assertions with a new suite of hydrodynamical 
simulations in scenarios where galaxy formation at high 
redshift preheats the intergalactic gas before it collapses 
to form clusters ( Bialck, Evrard & Mohr 2000 ; Mohr et al. 
in prep). However, most importantly, we emphasize that 
because of the sensitivity of X-ray surveys to evolution, we 
have only used those clusters which produce enough pho- 
tons to measure an emission weighted mean temperature. 
In this case, one can directly extract the minimum tem- 
perature Tum{z) of detected clusters as a function of red- 
shift. Correctly interpreting such a survey requires map- 
ping Tum{z) Miim{z) using the mass-temperature rela- 
tion; the evolution of the mass-temperature relation is less 
sensitive to the details of preheating than the luminosity- 
temperature relation. Thus, in a survey constructed in this 
manner, it should be possible to disentangle the cosmolog- 
ical effects from those caused by the evolution of cluster 
structure. 

Cluster Mass Function. In our treatment, we have re- 
lied on the mass function inferred from large scale numer- 
ical simulations of Jenkins et al. (2000). Although we do 
not expect the results presented here to change qualita- 
tively, changes in dN/dM by upto the quoted accuracy of 
~ 30% could affect the exact shape of the likelihood con- 
tours shown in Figures H and ^. It is important to test the 
scaling of the mass function with cosmological parameters 
in future simulations. We have further ignored the effects 
of galaxy formation and feedback on the limiting mass. In 
principle, the relation between the cluster SZE decrement 
and virial mass in the lowest mass clusters could be af- 
fected by these processes. In addition, the dependence of 
both the SZE decrement and the X-ray flux likely exhibits 
a non-negligible intrinsic scatter. The SZE decrement to 
virial mass relation is found to have a small scatter in 
numerical simulations (Metzler 1998), and to cause a neg- 
ligible increase in the total cluster yields (Holder et al. 
1999). However, the presence of scatter could effectively 
lower the limiting masses in our treatment of the X-ray 
survey. 

Local Cluster Abundance. Perhaps the most critical as- 
sumption is that the local cluster abundance is known to 
high accuracy. We have used this assumption to determine 
CTg, i.e. to eliminate one free parameter - effectively assign- 
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Fig. 12. — The middle panels show the likelihood contours for a fixed h = 0.65 in the X— ray survey, as in Figure ^ The upper and lower 
panels show the deviations in the contours caused by either a ±5% or a ±10% offset in the M — T normalization. 
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Fig. 13. — The middle panels show the likelihood contours for a fixed h = 0.65 in the X— ray survey, as in Figure^ The other two panels 
show the deviations caused by an offset in the slope M oc T'* . 



ing "infinite weight" to the cluster abundance near z = 0. 
This approach is appropriate for several reasons. The cos- 
mological parameters make little difference to the cluster 
abundance at z w 0, other than the volume being propor- 
tional to h~^. Similarly, the study of local cluster masses 
is cosmologically independent (upto a factor of h). In a 
10** square degree survey, we find that the total number of 
clusters between < z < 0.1, down to a limiting mass of 



2 X 10^*h~^ Mq is « 2500; with a random error of only 
±2%. We have experimented with our models, assuming 
that the normalization at z = is incorrectly determined 
by a fraction of 2%. In Figure |lj, we show the shift in the 
usual likelihood contour in the X-ray survey, caused by 
errors in the local abundance at this level. As the figure 
shows, the shift is relatively small (by about the width of 
the la region). In similar calculations with errors of ±4%, 
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we find shifts that are approximately twice as significant. 
We conclude that for our normalization procedure to be 
valid, the local cluster abundance has to be known to an 
accuracy of about 10%. 

Although such an accuracy can be achieved by only 
600 nearby clusters (which can be provided, for ex- 
ample, by an analysis of the SDSS data or perhaps the 
2MASS survey), it is interesting to consider a different ap- 
proach, where ug is treated as another free parameter in 
addition to flm, and w. The result of such a calcula- 
tion over a 4-dimensional grid is displayed in Figure 
This figure shows the likelihood contours along the slice 
h — 0.65 through this parameter space, but in projection 
along the erg axis: to be compared directly with the middle 
panel of Figure ^. Allowing as to vary results in a range 
of values 0.70 < cts < 0.97, and considerably expands the 
allowed likelihood region. The shape of the contours stay 
nearly unchanged, but their widths along the f2„j direc- 
tion expand by approximately a factor of 4, and their 
lengths along the w direction increase by about a factor of 
2. We conclude that our constraints would be significantly 
weakened without the local normalization (but would still 
be potentially useful when combined with other data; see 
below) . 

More General Cosmologies. In section |^, we restricted 
our range of models to flat CDM models. We find that 
the redshift distribution of clusters in open CDM models 
typically resembles that in models with high w. This is 
demonstrated in Figure ^: both in the w — —0.2 and the 
OCDM model, the redshift distributions are flatter and 
extend to higher z than in ACDM. We find that OCDM 
models with suitably adjusted values of flm and h are typ- 
ically difficult to distinguish from those with w — 0.5, 
but the flat shape of dN/dzdVl makes OCDM easily dis- 
tinguishable from ACDM. Note that open CDM models 
appear inconsistent with the recent CMB anisotrop y data 



conditions could reveal new degeneracies, and will be stud- 
ied elsewhere. 



from the Bo o merang and Maxima experiment s (e.g. Langc 
et al j 2000| , [White, Scott fc Pierpaoh 2000| , [Bond ct al 
2000^ A broader study of different cosmological mod- 
els, including those with both dark energy and curvature, 
time-dependent w, and those with non-Gaussian initial 
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Fig. 15. — Likelihppd contours for a fixed h = 0.65 in the X— ray 
survey, as in Figure Bi however, we here considered erg as a free pa- 
rameter, rather than fixing its value based on the local abundance. 
For the range of H and w shown, models with a likelihood better 
than 3(T take on values between 0.7 < erg < 0.97. 



6.3. Clusters versus CMB Anisotropy and High-z SNe 

A useful generic feature of the likelihood contours pre- 
sented here is their difference from those expected in CMB 
anisotropy or Supernovae data. Two different cosmolo- 
gies produce the same location (spherical harmonic index 
^pcak) for the first Doppler peak for the CMB tempera- 
ture anisotropy, provided they have the same comoving 
distance to the s urface of last scattering (cf. Wang & 



Steinhardt 199^ , [White 1998| , [Huey et al. 19991) . Note 



that this is only the most prominent constraint that can 
be obtained from the CMB data, with considerable more 
information once the location and height of the second and 
higher Doppler peaks are measured. Similarly; the appar- 
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Fig. 14. — The middle panel shows the likelihood contours for a fixed h = 0.65 in the X— ray survey, as in Figure joj The left and right 
panels show the deviations in the contours caused by a ±2% offset in the local cluster abundance determination. 
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ent magnitudes of the observed SNe constrain the luminos- 
ity di stance dr,(z) to < z ^ 1 (Schmidt ct aL 1998, Perl- 
mutter et aL 1999). In general, both of these types of 
observations will determine a combination of cosmological 
parameters that is different from the cluster constraints 
derived here. 
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Fig. 16. — Likelihood contours for a fixed h = 0.65 as in Fig- 
ure pi but zooming in for clarity. Also shown are combinations of 
w and Qrn that keep the spherical harmonic index £ of the first 
Doppler peak in the CMB anisotropy data constant to within ±1% 
(dashed lines); and combinations that keep the luminosity distance 
to redshift z = 1 constant to the same accuracy. 

In Figure |l^, we zoom in on the relevant region of the 
flrn — w plane in the X-ray survey, and compare the clus- 
ter constraints to those expected from CMB anisotropy 
or high-z SNe. The three dashed curves correspond to 
the CMB constraints: the middle curve shows a combina- 
tion of fljn and w that produces the constant ^pcak ~ 243 
obtained in our fiducial ACDM model (using the fitting 
formulae from White 1998 for the physical scale fcpcak); 
the other two dotted curves bracket a ±1% range around 
this value. Similarly, the dotted curves correspond to the 
constraints from SNe. The middle curve shows a line of 
constant at z = 1 that agrees with the ACDM model; 
the two other curves produce a c?l that differs from the 
fiducial value by ±1%. As the figures show, the lines of 
CMB and SNe parameter degeneracies run somewhat un- 
favorably parallel to each other; however, both of those de- 
generacies are much more complementary to the direction 
of the parameter degeneracy in cluster abundance stud- 
ies. In particular, the maximum allowed value of w, using 
both the CMB or SNe data, is w « —0.8; while this is 
reduced to w ^ —0.95 when the cluster constraints are 
added. Note that in Figure |l^, we have assumed a fixed 
value of h = 0.65; however, we find that relaxing this 
assumption does not significantly change the above con- 
clusion. The CMB and SNe constraints depend more sen- 
sitively on h than the cluster constraints do: as a result, 
the confidence regions do not overlap significantly even in 
the three-dimensional (w, fim, h) space. 

The high complementarity of the cluster constraint to 
those from the other two methods can be understood based 



on the discussions in S 4.1. To remain consistent with the 



CMB and SNe la constraints, an increase in w must be 
coupled with a decrease in flm] however, both increasing 
w and lowering ri,„ raises the number of detected clusters. 
To keep the total number of clusters constant, an increase 
in w must be balanced by an increase in fl^- Note that 
this statement is true both for the SZE and the X-ray 
surveys. Combining the cluster constraints with the CMB 
and SNe la constraints will therefore likely result in im- 
proved estimates of the cosmological parameters, and we 
do not expect this conclusion to rely on the details of the 
two surveys considered here. 

7. CONCLUSIONS 

We studied the expected evolution of galaxy cluster 
abundance from ^ z < 3 in different cosmologies, in- 
cluding the effects of variations in the cosmic equation of 
state parameter w = p/p. By considering a range of cos- 
mological models, we quantified the accuracy to which firm 
w, and h can be determined in the future, using a 12 deg^ 
Sunyaev-Zel'dovich Effect survey and a deep 10** deg^ X- 
ray survey. In our analysis, we have assumed that the 
local cluster abundance is known accurately: we find that 
in practice, an accuracy of ^ 5% is sufhcient for our results 
to be valid. 

We find that raising w significantly flattens the redshift- 
distribution, which can not be mimicked by variations in 
either flm, h, which affect essentially only the normaliza- 
tion of the redshift distribution. As a result, both surveys 
will be able to improve present constraints on w. In the 
— w plane, both the SZE and X-ray surveys yield con- 
straints that are highly complementary to those obtained 
from the CMB anisotropy and high-z SNe. Note that the 
SZE and X-ray surveys are themselves somewhat comple- 
mentary. In combination with these data, the SZE survey 
can determine both w and U,m to an accuracy of « 10% 
at 3<T significance. Further improvements will be possi- 
ble from the X-ray survey. The large number of clusters 
further alleviates the degeneracy between w and both f2m 
and h, and, as a result, the X-ray sample can determine w 
to « 10% and to w 5% accuracy, in combination with 
either the CMB or the SN data. 

Our work focuses primarily on the statistics of cluster 
surveys. We have provided an estimate of the scale of 
various systematic uncertainties. Further work is needed 
to clarify the role of these uncertainties, arising especially 
from the analytic estimates of the scaling of the mass lim- 
its with cosmology, the dependence of the cluster mass 
function on cosmology, and our neglect of issues such as 
galaxy formation in the lowest mass clusters. However, 
our findings suggest that, in a flat universe, the cluster 
data lead to tight constraints on a combination of and 
w, especially valuable because of their high complementar- 
ity to those obtained from the CMB anisotropy or Hubble 
diagrams using SNe as standard candles. 
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